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Summary 


This paper contains a collection of some results of four individual studies presenting 
calculated numerical values for airfoil aerodynamic stability derivatives in unseparated 
inviscid incompressible flow due separately to angle-of-attack, pitch rate, flap deflection, 
and airfoil camber using a discrete vortex method. Both steady conditions and oscillatory 
motion were considered. Variables include the number of vortices representing the airfoil, 
the pitch axis / moment center chordwise location, flap chord to airfoil chord ratio, and 
circular or parabolic arc camber. Comparisons with some experimental and other 
theoretical information are included. The calculated aerodynamic numerical results 
obtained using a limited number of vortices provided in each study compared favorably with 
thin airfoil theory predictions. Of particular interest are those aerodynamic results calculated 
herein (such as induced drag) that are not readily available elsewhere. 


Introduction 

Calculated numerical values for aerodynamic coefficients and stability derivatives 
for some airfoils in unseparated inviscid incompressible flow using a discrete vortex method 
consisting of a limited number of two-dimensional vortices were obtained for both steady 
conditions and sinusoidal oscillatory motion. Derivatives were calculated individually with 
respect to angle-of-attack, pitch rate, and flap deflection. In addition the effect of airfoil 
camber in steady flow was examined separately. The aerodynamic derivatives due to each 
of the four variables were examined in considerable detail and the information is presented 
here as a collection of the results of the four individual studies labeled as Section A through 
D in the order in which the effort was undertaken. The collection is believed to address 
some of the aerodynamic interest in airfoils at low speed and oscillatory frequencies 
relevant to vehicle dynamic stability. 

A number of theoretical studies of airfoils in steady and oscillatory conditions exist in 
the published literature over the past 75 years (for oscillatory examples see textbooks like 
reference 1, numerous NACA/NASA Reports, AIAA Journal Articles, AGARD publications 
like reference 2, and many others). These studies examine the problem for a range of 
conditions including subsonic and supersonic Mach Numbers as well as for incompressible 
flow. Nearly all mathematical developments were restricted to the linear range of small 
angles and rates and address the airfoil flow condition at the trailing edge (the Kutta 
condition). The discrete vortex method circumvents the Kutta condition by substituting a 
boundary condition at the three-quarter-chord location that yields equivalent aerodynamic 
results. The purpose of this effort, which was part of a larger exploratory investigation 
including both airfoils and wings, was to apply the discrete vortex method to assess the use 
of a limited number of vortices to yield realistic values of the aerodynamic derivatives. 



Symbols 


A amplitude 

c airfoil chord, ft. 

c f flap chord, ft. 

c d section drag coefficient, d / .5 pV~ c 

c h section flap hinge-moment coefficient, h / . 5pV 2 [c f j 

c e flap lift coefficient, £l.5pV~c 

c, section lift coefficient, l/.5pV~c 

c m section pitching-moment coefficient, m/.5pV 2 c 2 

c p pressure coefficient, A/?/. 5 pV~ 

d airfoil drag force, lbs. per unit span 
h flap hinge-moment, ft. -lbs. per unit span 
k reduced frequency parameter, coc/2V 

K n individual vortex circulation value, ft 2 / sec. 
t flap lift, lbs. per unit span 

l airfoil lift force, lbs. per unit span 

m airfoil section pitching moment, ft. -lbs. per unit span 
N number of vortices used 

p mass density of air, slugs / ft 3 

V velocity, ft. /sec 

t time, sec. 

a y e 

Z camber parameter, — (circle) or — (parabola) 

c c 

a airfoil angle of attack, radians 

C 1 airfoil pitch-rate, radians per sec. 

d flap deflection, radians 

CCf flap angle of attack, radians 

q f flap pitch-rate about flap leading edge, radians per sec. 
CO circular angular velocity, radians/sec. 

6 pitch angle, radians 

O indicial function 

0 slope of the mean line, radians 

(j) e trailing edge angle, radians 


2 



deb 

c b = — y ~ 

a da 


cb = 

a 


d£b_ 

ac 
a 

2V 


c b = 

q 


d c b 
d V- 

2V 


dc b ri _ dc b ... d£b 

’af~ daf b af~ qf~ d & 

2V 2V 



dcb 

dd 


Oscillating Conditions: 

cb = cb -k 2 cb cb = cb + cb 

S of qf 8 qf «/ 

where b = l. , £, h, m 


Subscripts: 

0 amplitude of sinusoidal motion 

a due to angle of attack 

q due to pitch rate 

ss steady state 
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Theory 

A thin symmetrical airfoil is chosen to be represented by its chord that is divided into 
a number of segments. A vortex is located at the one-quarter position of each segment and 
a control point at the three-quarter position. The velocity generated by the sum of all 
vortices at each control point is set equal to the external wind velocity normal to the airfoil 
surface. Solution of the set of simultaneous equations yields values for the unknown 
circulations and thus the lift for steady conditions. For time-dependent calculations use is 
made of starting vortices positioned at a quarter-chord aft of the trailing edge at time equals 
zero. As time increases vortices are shed downstream in the wake. Calculation of the lift on 
the airfoil surface at each new wake location yields the response to a step-input. From this 
indicial function, additional calculations can provide the frequency response from which the 
aerodynamic derivatives as a function of reduced frequency, k, can be obtained. 
Calculations are made separately for angle-of-attack and pitch-rate, and when summed the 
result corresponds to wind-tunnel sinusoidal oscillatory measurements (see Appendix A). 


3 



Presentation Of Results 


The four individual studies are identified briefly as follows: 

Section A - Steady aerodynamic derivatives for angle of attack and pitch rate 
Section B - Oscillatory aerodynamic derivatives for angle-of-attack and pitch rate 
Section C - Steady and oscillatory aerodynamic derivatives due to flap deflection 
Section D -Steady aerodynamics due to circular and parabolic arc camber 


Results and Discussion 

The calculated aerodynamic results presented in each section compared favorably 
with thin airfoil theory predictions. Of particular interest are those aerodynamic results 
calculated herein that are not readily available elsewhere. These are worth noting and 
include: (1) the load distribution across the chord for steady pitch rate, (2) values for the 
induced drag coefficient for oscillatory motion, (3) airfoil and flap motions that result in lift 
and pitching moment versus displacement traces having unusual shapes, and (4) slight 
differences in the aerodynamics between circular and parabolic arc camber. 


Section A - Steady aerodynamic derivatives for angle of attack and pitch rate 

Consider the airfoil section to be symmetrical and thus represented by its chord. 
Calculated values for the angle-of-attack and pitch-rate derivatives for different numbers of 
vortices equally spaced across the chord are tabulated respectively in Tables 1 and 2 for 
five different pitch-axis/moment-center locations on the chord. The values are exact and 
presented in the form of fractions as calculated. Thin airfoil theory values (first given in 
1928 by Glauert in reference 3 and converted here to NASA nomenclature) are also 
included for comparison and plotted in figure 1 to illustrate the parabolic curve for pitch 
damping derivative c m . Calculated values herein for Ci , c m , and C; are independent 

rn q L a fn a L q 

of the number of vortices employed and identical to those of thin airfoil theory on figure 1 . 
Values for c m however approach thin airfoil values when a reasonable number of vortices 

rn q 

are used to represent the airfoil as illustrated by figures 2 and 3. Figures 4, 5, and 6 are 
provided to indicate the load distribution across the airfoil chord due to angle of attack and 
due to pitch rate. Of particular interest is the zero value of the moment derivative due to 
pitch rate at the mid-chord position resulting from the up-loading on the front as well as the 
rear portions of the airfoil. Finally, airfoil induced drag predicted by the interaction of the 
various vortices for both angle-of-attack and pitch-rate conditions was found to be zero as 
expected and is illustrated in the example on figure 7 (for additional comments on drag see 
Appendix B). 
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Table 1. - Aerodynamic Derivatives Due to Angle of Attack 



Moment Center Chordwise Location 

Leading 

Edge 

Quarter 

Chord 

Mid 

Chord 

Three-Quarter 

Chord 

Trailing 

Edge 

c l 

'a 

C m 
11 'a 

c l 

'a 

C m 

m a 

c l 

'a 

c m 

IM a 

c l 

'ex 

c m 

■■■a 

c l 

'ex 

C m 
11 'a 

Mathematical Exact Value 

Number of Vortices 

1 

2jt 

1 

“T" 

2jt 

0 

2jt 

1 

+ -r-3t 

2 

23T 

+1 ax 

23t 

3 

+ "T - 3t 

2 

2 

2jx 

1 

~T n 

23T 

0 

23t 

1 

+ 3T 

2 

23T 

+ 1 3T 

23t 

3 

+ K 

2 

3 

2k 

i 

~T n 

23T 

0 

23t 

1 

+ 3T 

2 

23T 

+ 1 3T 

23t 

3 

+ 3T 

2 

D 

2k 

i 

-y 31 

23t 

0 

23t 

1 

+ -r-3t 

2 

23T 

+ 1 3T 

23T 

3 

+ -T“3t 

2 

00 

* 

2jz 

i 

-y 31 

23T 

0 

23T 

1 

+ -r-3t 

2 

23T 

+ 1 3T 

23t 

3 

+ -T“3t 

2 

Computer Calculated 

1 Values 

8 

2k 

-.50003t 

23t 

0 

2k 

.50003T 

2k 

1 .OOOOjc 

2it 

EU2!K 


* Thin airfoil theory values 


Table 2. - Aerodynamic Derivatives Due to Pitch Rate 



Pitch-Axis / Moment-Center Chordwise Location 

Leading 

Edge 

Quarter 

Chord 

Mid 

Chord 

Three-Quarter 

Chord 

Trailing 

Edge 

C| q 

c m q 

C| q 

c m q 

C| q 

c m q 

C| q 

c m q 

C| q 

c m q 

Mathematical Exact 

: Value 

Number of Vortices 

1 

33t 

3 

~T n 

2jt 

0 

K 

1 

+ T * 1 

4 

0 

0 

— 3T 

3 

2 

3jt 

15 

~16 n 

2jt 

3 

~16 n 

K 

1 

+ 16 n 

0 

3 

~16 n 

-3T 

15 

“16^ 

3 

33t 

35 

“36™ 

23T 

8 

~36 n 

K 

1 

+ 36 3t 

0 

8 

~ 36 n 

— 3T 

35 

_ 36 31 

D 

33T 

63 

~FT n 

64 

23t 

15 

-FT n 

64 

K 

1 

+ FT n 

64 

0 

15 

- 77 3t 

64 

— 3T 

63 

-7T n 

64 

00 

* 

33t 

-K 

23T 

1 

K 

0 

0 

1 

— 3T 

-K 

Computer Calculatec 

1 Values 

8 

3jt 

-.9961 k 

2jt 

-.24613T 

K 

.0039k 

0 

-.246l3t 

-K 

-,996l3t 


* Thin airfoil theory values 
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Section B - Oscillatory aerodynamic derivatives for angle of attack and pitch rate 

Consider the airfoil to be symmetrical and thus represented by its chord. 

Calculations were made considering the airfoil to be in rectilinear motion at time equals 
zero and then a step input is applied in angle of attack and separately in pitch rate. A 
starting vortex is positioned at c/4 aft of the airfoil trailing edge at t=0 that moves 
downstream with increasing time. Figures 8 and 9 illustrate the vortex arrangements that 
were used in the calculations for the wake at various positions downstream. The resulting 
indicial function developed was then used to determine the frequency response. As in the 
case for steady motion, calculations were made for five pitch-axis locations equally spaced 
across the chord from the leading edge to the trailing edge. Figures 10 and 1 1 give the lift 
and pitching moment results for a step input in angle of attack as a function of downstream 
location. Figure 12 gives the corresponding induced drag function resulting from the 
interaction of the vortices in the wake with those on the airfoil. The equation on figure 12, 
although not an indicial function, can be used to estimate the induced drag coefficient with 
wake positions further downstream. An examination of the figure shows the maximum 
induced drag value occurs just before the wake has proceeded one chord length 
downstream. As time increases and the wake continues to move downstream the induced 
drag force decreases until the steady-state condition is reached at which time the induced 
drag is zero. 

Using the indicial functions developed from figures 10 and 11, calculations were 
made to obtain the corresponding lift and pitching-moment frequency response results. 
Figures 13 and 14 give the lift and pitching-moment derivatives Cj ,Cj. ,c m , and c m . as 

a function of reduced frequency parameter k. Figures 15 and 16 are included here to 
illustrate the use of figures 13 and 14 to obtain typical responses to sinusoidal plunging 
oscillations of a thin two-dimensional airfoil. These figures show lift-displacement and 
moment-displacement traces for a given value of reduced frequency. In the case of the 
induced drag, figure 12 shows E c i a is not an indicial function, however the square root of 
Erf a labeled ^d a in figure 17 can be considered as such. Frequency response results 
permit establishing a sinusoidal time history trace of the function. Squaring this result 

provides a time history of the induced drag parameter that is twice the motion frequency 
and from which the induced-drag — displacement trace can be constructed. Note that to 
obtain a drag comparison with experimental data, the profile drag contribution must be 
added to the induced drag component shown. 

Calculated results for the lift, pitching moment, and induced drag due to a step input 
in pitch rate were also determined. From this information, frequency response results for 
the lift and pitching moment derivatives shown on figures 18 and 19 were obtained as a 
function of reduced frequency parameter k. Step-response results obtained for the induced 
drag function due to pitch rate are given in figure 20 for completeness, however further 
effort with this induced drag information was not undertaken. 

Combination derivatives representative of low-speed wind-tunnel aerodynamic 
measurements of an airfoil pitching about a point on the airfoil chord were obtained by 
appropriately summing the angle-of-attack and pitch-rate derivatives. Figures 21 and 22 
provide these derivative values for the five different pitch-axis/moment-center locations on 
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the chord. An examination of the numerical results in general show a reduction in the 
derivative magnitude with an increase in oscillation frequency. 

Some comparisons of values from calculations made herein using 4 vortices to 
represent the airfoil with other theoretical values and experimental measurements are 
presented in figures 23 thru 27. Figures 23 and 24 give comparisons of the airfoil indicial-lift 
function with Wagner’s prediction and figure 24 with Theodorsen’s prediction. (See 
reference 5.) Comparison of the combined derivative values calculated herein with wind- 
tunnel measurements given by Halfman in reference 4 converted to the nomenclature of 
this paper and adjusted by the inverse of the experimental measured lift-curve-slope ratio 

c l 

— — of 0.813 for steady conditions is given in figures 25, 26 and 27. The experiment was 
2 it 

conducted on a 1-foot chord, 2-foot span NACA 0012 airfoil mounted between two large 
vertical fairings in the 5- by 7- foot test section of the Massachusetts Institute of Technology 
flutter tunnel. The supporting mechanism could provide sinusoidal oscillatory vertical 
motion (±linch and ±2 inches) and angular rotations about the 37% chord position either 
independently or in combination for a reduced frequency range for k from about .05 to 
about .45. Considering the differences existing between the mathematical models and the 
wind-tunnel experimental set-up, the data comparisons shown in the figures appear 
reasonable. 


Section C - Steady and oscillatory aerodynamic derivatives due to flap deflection 

The airfoil is considered symmetrical so that it can be represented by its chord. The 
chord is divided into a number of segments according to the number of vortices chosen to 
represent the airfoil. Vortices and control points were positioned on the chord line according 
to the one-quarter — three-quarter rule for each segment. The hinge line for flap deflection 
was placed at a segment leading edge. This arrangement corresponds to a plain flap with 
gap sealed. Calculations for selected values of the ratio of flap chord to airfoil chord were 
chosen to cover the range from 0 tol . All angle and rates are assumed small. Steady 
conditions are considered first followed by oscillatory motion. 


Calculated values for the airfoil aerodynamic derivatives and c m using different 




numbers of vortices representing the airfoil are given in figure 28. The ordinate on 


figure 28 also can be interpreted as 


Aa 

A<5 


2 it 


As a consequence the total lift coefficient with 


flap deflected can be written as C/ = C/ 


A a 

a + — o 

A <5 


since the method of superposition 


applies. Calculated values for steady conditions of lift derivatives for the airfoil and 
separately for the portion of the lift derivative carried on the flap due to angle of attack is 
given in figure 29a. Shown on figure 29b are the corresponding lift derivatives due only to 
flap deflection. Some interesting comparisons can be observed from an examination of the 

c f 1 

values in the two figures. For example, consider the derivative values for — = — . For the 

c 2 
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flap deflected (fig. 29b), 50% more lift is produced on the flap compared to the lift when the 
flap is in the undeflected condition (fig 29a). Lastly, deflecting the half-chord flap produces 
80% of the lift produced by deflecting the complete airfoil the same amount. Calculation of 
the hinge-moment derivatives, Cu and C/,„ were made using 8 vortices representing the 

n a n o 

airfoil and the results are given in figure 30. Derivative values were obtained for selected 
values of the ratio of flap chord to airfoil chord. The total hinge-moment coefficient can be 
written as 

c h ~ ( 'h (f Cl + ('h () <5 


Also shown on figure 30 for comparison purposes are values obtained from reference 5. 


For oscillatory motion, calculation of the aerodynamic response due to a step input 
in flap deflection (indicial response) is similar to that used to obtain the response for the 
airfoil in Section B. A starting vortex is located at a distance c/4 aft of the trailing edge 
and consists of all vortices used to define the airfoil with deflected flap. As the wake moves 
downstream with increasing time, additional vortices are inserted on the airfoil and in the 
wake to satisfy the boundary conditions. Aerodynamic calculations are carried out 

separately for a step input in flap attitude (flap angle of attack, CLf) and separately in flap 


angular velocity (flap pitch-rate, qj). Figures 31, 32, 33, 34, 35, 36, and 37 are presented 

Vt _ 

giving some aerodynamic derivatives for the starting condition = U and the steady 

cl 2 


state condition 


Vt 


c/2 

dependency involved. 


= 00 . The difference in the curves is an indication of the time 


Once the indicial functions are determined the frequency response can be obtained. 
Figures 38 and 39 give the lift derivatives as a function of reduced frequency parameter k 
for two flap sizes. The figures give a comparison of the flap data with that of the complete 
airfoil both oscillating at the same frequency. It is interesting to note that the derivative 
values for the half-chord flap oscillating are about three-quarters of those for the complete 
airfoil oscillating. 


Since the method of superposition holds for combining the lift and moment 
aerodynamic derivatives due individually to the CL andg motions to yield values for 
sinusoidal oscillations, superposition can be applied to different combinations of amplitude 
and frequency of airfoil and flap motion. Figures 40, 41 , 42, and 43 give several example 
lift-displacement traces that were calculated. Reference 6 presents a number of 
comparisons of airfoil and flap data from wind tunnel measurements with a different 
analytical computation showing loops like those on figures 41 and 42. 


Section D - Steady aerodynamics for circular and parabolic arc camber 

A limited number of calculations using a discrete vortex method in steady inviscid 
incompressible flow were made to obtain the effect of camber on the lift and pitching 
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moment of a two dimensional airfoil represented by its mean-line. Calculations were made 
for circular arc camber and also for parabolic arc camber. The latter camber was 
undertaken since NACA Four-Digit Airfoils use mean lines consisting of parabolic arcs thus 
permitting comparison of calculated and experimental results. Maximum camber values of 
1 , 2, and 4 percent of the chord along with maximum camber locations of 40, 50, and 60 
percent of the chord were considered. The geometry involved and calculated values of lift 
and pitching moment are given and discussed. To permit calculation for different chordwise 
locations of maximum camber, the airfoil mean line is divided into two parts at the point of 
maximum camber. Each part is then divided into a number of segments with a vortex and 
control point assigned to each segment. For each segment, the one-quarter — three-quarter 
rule is applied. Since displacement of the mean line from the chord line is small for small 
camber, the vortices and control points are placed on the airfoil chord instead of the curved 
mean line to simplify the geometry involved. However, the external velocity at each control 
point to be canceled by the vortex system is that at the mean line since that is the local 
angle of attack of the surface at that point. The external velocity at each control point can 
be expressed in terms of the trailing edge angle (f) e . Solution of the matrix equation yields 
values of the unknown circulations as a function of the trailing edge angle <p e . Using the 
relation between trailing edge angle (f) e and camber parameter z along with the Kutta- 
Joukowski equation gives the lift and pitching moment coefficients as a function of z where 
Z is the maximum displacement of the mean line from the chord line in fractions of the 
chord. Note that for these calculations, the airfoil chord line must remain aligned with the 
relative wind (i.e., airfoil angle of attack remains zero). Three factors of interest in the use 
of camber are: 

1) the amount of camber 

2) the chordwise location of the maximum camber point 

3) the maximum airfoil thickness 

Other influences such as the equation of the camber line, Reynolds Number, Mach 
Number, etc., are also of concern but the three main factors were important enough that 
the NACA Four-Digit Airfoil Series identifies them in the airfoil designation. For example the 
NACA 2412 section specifies the airfoil has 2% camber (i.e., the maximum displacement of 
the mean line from the chord line expressed in fractions of the chord), the maximum 
displacement is at 40% of the chord from the leading edge, and the 1 2 represents the 
maximum thickness in percent of the chord (for additional background discussion on 
Camber see Appendix C). 

Calculations using the geometry given in figure 44 were made to obtain the lift and 
pitching moment coefficients due to circular arc camber with the maximum camber point at 
the 50% chord location. Calculations were made using 2 and 4 vortices to represent the 
airfoil. Table 3 lists the calculated aerodynamic coefficients in terms of the trailing edge 
angle (j) e and also in terms of the camber parameter z . 
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Table 3 - Aerodynamic Coefficients For Max-Camber At Mid- Chord 


Number 

Of 

Vortices 

Max 

Camber 

Location 

c l 

c ' m c/2 

C}n c / 4 

2 

,50 c 


m Pe 

i 

ON OJ 
1- 

4 

,50 c 

Jt(p e 








2 

,50 c 

4 JTZ 


~4 Jrz 

4 

,50 c 

4 jxz 

16® 

16 - 

00 

,50 c 

4 jiz 

0 

-nz 


Note that the table indicates that the value of lift coefficient is independent of the number of 
vortices used to represent the airfoil.. The value of the pitching moment coefficient, 
however, depends on the number of vortices representing the airfoil. The table also lists the 
limiting values. The zero value of the moment at the mid-chord location indicates this is the 
position of the center of pressure. McCormick in Reference 7 presents results using a large 
number of vortices illustrating the limiting value about the mid-chord is zero. The procedure 
used for the calculations herein is similar to that of Reference 7 with similar results. 
Differences, however, exist in the development of the camber geometry presented in figure 
44c. 


To illustrate the influence of altering the chordwise location of maximum camber, 
some additional calculations were made for maximum camber at the 40% and 60% chord 
location. For each condition, two circular arcs of different radius were chosen such that they 
were tangent at the point of maximum camber and that the airfoil chord was aligned with 
the relative wind. Figure 45 and 46 give the geometry used. Calculations were made using 
2, 4, and 8 vortices representing the airfoil and the results are given in Table 4. 
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Table 4 - Aerodynamic Coefficients For Max-Camber At Different Chord Positions 


Number 

of 

Vortices 

Maximum 

Camber 

Location 

Cl 

Cm C/4 

2 

.36c 

.8771 (4itz) 

-.5658(jtz) 

2 

,40c 

.9091 (4itz) 

-.6236(itz) 

2 

,44c 

.9434 (4itz) 

-.6651 (jtz) 

2 

,50c 

1 .0000 (4itz) 

-.7500(itz) 

2 

,56c 

1 .0638 (4utz) 

-.8457(itz) 

2 

,60c 

1.1111 (4jiz) 

-.9167(jtz) 





4 

,40c 

.9068 (4itz) 

-.7888(jtz) 

4 

,50c 

1 .0000 (4itz) 

-.9375(:tz) 

4 

,60c 

1.1236 (4jtz) 

-1.1 221 (jtz) 





8 

,40c 

.9063 (4itz) 

-.8316(itz) 

8 

,50c 

1 .0000 (4jiz) 

-.9837(jtz) 

8 

,60c 

1.1286 (4jiz) 

-1.1744(jtz) 





00 

,40c 

.9063 (4itz) 

-.8479(jtz)* 

00 

.50c 

1 .0000 (4jiz) 

-I.OOOO(jtz) 

00 

,60c 

1.1286 (4itz) 

-1 .1907(jtz)* 






* Estimated using increment at mid chord 


For these calculations the numerical values were expressed in digital form using eight 
decimal places. An examination of the tabulated values indicates both the effect of camber 
location and the number of vortices used. A graphic illustration is given in figure 47 and 
sketches of the chordwise load distributions are given in figure 48. Comparison of the 

ordinate numerical values of figure 48 cannot be made because the value of (p e is different 

for the three distributions, however the shift in loading is apparent with maximum camber 
location. An examination of figure 47 indicates that shifting the maximum camber location 
forward of the mid-chord decreases the lift and adds a pitch-up increment to the moment 
for a given vortex representation. Shifting the maximum camber point rearward of the mid- 
chord gives the opposite aerodynamic result for a given vortex representation. As indicated 
previously in the text, the lift coefficient was found to be independent of the number of 
vortices employed. This fact is also shown on figure 47. Of more significance however, was 
the large effect on the pitching moment of increasing the number of vortices used for the 
calculations. Increasing the number of vortices representing the airfoil added an increment 
to the pitch down moment. The values for the 50% maximum camber location indicated 
using 8 vortices provided a numerical value of the pitching moment coefficient within 2% of 
that using an infinite number of vortices. 

The camber geometry involved for NACA Four-Digit Airfoils is illustrated in figure 49 
and consists of two parabolic curves, one curve on each side of the maximum camber 
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position. (See reference 8.) Recall that a parabolic curve is the locus of points equidistant 
from a line called the directrix and a point called the focus. The equation for the curve was 
obtained from reference 10. Note that the origin of the axis system is midway between the 
directrix and the focus. The discrete vortex method places the vortices and control points 
representing the airfoil on the chord to simplify the geometry involved, however the flow 
velocity to be cancelled by the vortex system is that on the curved surface. Calculations 
were made for 2, 4, and 8 vortices for maximum camber at the mid chord. The numerical 
results obtained were the same as those obtained for circular arc camber. In addition 
calculations were made for maximum camber at 40% and 60% chord locations. Since 
reference 9 provides experimental measurements for maximum camber at the 40% chord 
location, a comparison of computed values for 2% and 4% camber with values from 
reference 9 is given in the following Tables. 


Table 5 - Comparison Of Calculated Values And Experimental Measurements 


Vortex Calculated 2%Camber 


Number Of Vortices 

C l 

Cm c/4 

8 

0.2278 

-0.05225 

00 

0.2278 

-0.05328 


Reference 9 Experimental Values 


Airfoil Section 

ci 

Cm d 4 

NACA 2408 

0.195 

-0.038 

NACA 2410 

0.20 

-0.050 

NACA 2412 

0.23 

-0.035 


Vortex Calculated 4% Camber 


Number Of Vortices 

Cl 

Cm cl 4 

8 

0.4556 

-0.10450 

00 

0.4556 

-0.10655 


Reference 9 Experimental Values 


Airfoil Section 

ci 

Cm d 4 

NACA 4412 

0.39 

-0.095 


The comparison of calculated and experimental values for the 2% camber information 
seem reasonable for the thickness ratios presented considering the difficulty of establishing 
the numerical values from the figures of reference 9. For the 4% camber information it 
appears that the calculated values overestimate the magnitude of the lift and pitch-down 
moment. 

The original camber derivation developed previously gave identical expressions for 
the aerodynamic lift and pitching moment coefficients for circular arc and parabolic arc 
camber. Since exact solutions would be expected to show some difference, the camber 
derivations were re-examined and an alternate procedure was developed for circular arc 
camber. Note that the origin of the axis system for the parabolic camber derivation on figure 
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49 was chosen to be at the vertex of the parabola on the Y-axis. For a similar situation to 
exist for circular arc camber the origin of the axis system must be moved from the center of 
the circle to a position a radial distance along the positive Y-axis. Using this set-up different 

relationships can be developed for parameter z and the trailing-edge angle (p e as shown 

in Appendix D. Using this new equation for circular arc camber, the new equations for lift 
and pitching moment coefficients are: 

c, =4jtz(l + 4z 2 } c mc/4 =-jzz{l + ^z 2 ) 

The factor ^1 + 4 z~ j indicates that for 2% camber, the increase in the aerodynamic 

coefficients is of the order of one-sixth of one percent of the original value. This small 
amount is the difference in the aerodynamic results between airfoils using circular arc and 
parabolic arc camber. 


Concluding Remarks 

Aerodynamic coefficients and stability derivatives for airfoils in unseparated inviscid 
incompressible flow were calculated using a discrete vortex method involving a limited 
number of two-dimensional vortices for both steady conditions and sinusoidal oscillatory 
motions at low speed and low frequencies. Comparison of calculated values for both steady 
state and oscillatory conditions were in fairly good agreement with those of thin airfoil 
theory. Limited comparisons of calculated numerical values herein with available wind- 
tunnel experimental measurements were considered reasonable. For oscillatory conditions 
induced drag values were calculated due to the interaction of the shed vortices in the wake 
with the vortices on the airfoil surface. For the steady condition induced drag is, of course, 
zero. Values of steady-state derivatives due to flap deflection were found to be in good 
agreement with previous published results. Aerodynamic lift for some different 
combinations of sinusoidal oscillatory motion of the airfoil and flap produced loops in the lift- 
displacement traces similar to those observed experimentally. Finally in steady flow, only 
very small differences were found to exist in the calculated aerodynamic results between 
airfoil mean-lines having circular or parabolic arc camber for small camber values. The 
position of the maximum camber point on the chord was found to have a significance 
influence on the values of the lift and pitching-moment coefficients. Comparison of 
calculated lift and moment magnitudes were found to overestimate somewhat the few 
experimental wind-tunnel measurements available. 
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Appendix A - Wind-Tunnel Measurements 


Consider a thin two-dimensional airfoil performing sinusoidal oscillations in pitch 
about an axis on the airfoil chord. Assume the air-stream is horizontal and in a direction 
parallel to the wind-tunnel centerline. If 0 designates the pitch angle of the airfoil relative to 
the wind-tunnel centerline, then the following expressions apply: 


9 = 9 0 sin rot 
9 = d 0 0 ) cos cot 
6 = - 9 0 u> 2 sin cot 


where 9 is the instantaneous pitch angle of the airfoil at any time t and 9 0 is the amplitude 
of the motion. The aerodynamic lift coefficient due to the oscillatory motion can be written 
as: 


cue 


qc 


q 


C|-C| a « +C la — + Q q — + C,A -(-) 


2V 


2V 


4 V 


Recognizing that 

a = 0 q = 9 

a = 9 q = 9 

Then upon proper substitution gives 

OJC (JOC coc 

C| = C| a 0 o sin rot + C| d 9 0 — cos rot + C| q 6 0 — cos rot - C\ q 0 O ( — ) 2 sin rnt 


Collecting similar terms and letting k = — gives: 

C| = ( C| a - k 2 C\ q ) 0 o sin cot + ( C| q + C| d ) k 9 0 cos cot 

Each term of this equation is a function of the motion amplitude and each is a function of 
reduced frequency. The first term is in-phase with the airfoil angular displacement 0 and the 
second term is out-of-phase with 0. An analogous development can be made for the 
pitching moment yielding: 

Cm = (Cm a ' k 2 Cm^r ) 0 o sin cot + (Cm^ + C m ^ ) k0 o COS rot 

An equation for drag coefficient similar to the lift and moment equations given above cannot 
be written since the induced drag is a nonlinear function of the a and q variables. 
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Appendix B - Comment On Airfoil Drag At Steady Flow Conditions 


It is interesting to observe that Jean Le Rond d’Alembert (a famous French 
mathematician born in 1717 and died in 1783) examined the airfoil drag problem for steady 
unseparated inviscid incompressible (potential) flow. The inviscid flow assumption specified 
that the friction drag component is zero. Potential flow specified that the flow leaves 
smoothly from the trailing edge along a streamline (the Kutta condition). Thus no separation 
and no discontinuities exist at the trailing edge. As a consequence, full pressure recovery 
occurs on the rear of the airfoil and thus pressure drag is zero. Since profile drag is the sum 
of the friction drag and pressure drag components, it is also zero. The remaining drag 
component that can exist is that of induced drag and it also is found to be zero. Thus 
occurs d’Alembert’s paradox — the two dimensional airfoil at an angle of attack produces a 
lift force but pays no penalty since drag is zero which is a quite unusual and unexpected 
result. d’Alembert worked to obtain a drag force rather than this unusual result for a 
number of years without success. Von Karman in reference 1 1 indicates zero drag is the 
correct result and it is the result given by the discrete vortex method. 


Appendix C - Background Discussion on Airfoil Camber 

Low speed aerodynamic characteristics for wings and airfoils are divided into 
planform effects and section effects. Planform effects involve variations in shape, aspect 
ratio, taper ratio, sweep angle, etc. The first useful developed wing theory was Prandtl’s 
lifting-line theory followed much later by vortex lattice and lifting surface schemes among 
others. Section effects were addressed by Max Munk (references 12 and 13) and others in 
the early 1920’s and resulted in the well known “thin airfoil theory”. In this theory the airfoil 
section is represented by its mean line that is the line from the leading edge to the trailing 
edge halfway between the upper and lower airfoil surfaces. Lift and pitching moment are 
associated with the shape of this line. Thickness distribution about the mean line provides 
the section shape and is responsible along with skin friction for the airfoil drag force. If the 
airfoil has a symmetrical section then the mean line is the chord line. For this situation the 
aerodynamic center is located at the one-quarter chord location, the lift curve slope is 2 jt, 
and the pitching moment coefficient is zero. For the mean line displaced from the chord line 
both lift and pitching moment are obtained. Curvature of the mean line is usually 
designated by an equation and denoted by the maximum vertical displacement from the 
chord and its chordwise station aft of the leading edge. For airfoil section lift and pitching 
moment, angle of attack and curvature effects are additive for steady conditions. Because 
of the presence of the boundary layer, some difference exists between theory and 
experiment. Section drag is usually obtained experimentally. Over the years a large number 
of airfoil sections have been designed and tested experimentally in wind tunnels. A 
collection of some of this information is given in a book entitled Theory of Wing Sections by 
Abbott and Von Doenhoff (ref. 9). 


Appendix D - Alternate Circular-Arc Camber Derivation (max camber at 50% chord) 

The origin of the axis system for the parabolic camber derivation (figure 49) was 
chosen to be at the vertex of the parabola on the Y-axis. For a similar situation to exist for 
circular arc camber, the origin of the axis system must be moved from the center of the 
circle to a position a radial distance along the positive Y-axis. Using this set-up, 
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relationships can be developed between the camber parameter z and the trailing-edge 
angle (j) e . The equation for a circle with the axis system origin at the center of the circle is: 

2,2 n 2 

x + y = R 

The equation of a circle with the axis-system origin displaced vertically with the X-axis 
tangent to the circle is: 

x 2 + (y + R) 2 = R 2 
x 2 + y 2 + 2 Ry' + R 2 = R 2 

2 l2 1 } I 

x + y = -2 Ry 


Taking the derivative 


2 xdx + 2 y'dy' = -2 Rdy 
xdx + y'dy = -Rdy' 
xdx = -(y + R)dy' 
dy' x 
dx y + R 


Replacing R gives 


x 2 + y' 2 \ -2y' z + x z + y' z x z -y 


2 . ..2 . ,il 2 .2 


Now dx ~ * an ^ ~ ^ f° r small values of (f> then 


x x , , c a y? 

When q> = cp e , y =y e , x = x e = — and z = — = — then 

2 c c 


9 Je A Je 

c/2 _ C 

/ i \ 2 / t \2 

l-i 1-4 — 

\c/2 \ c 
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4 z 

1-4 z 2 


0 ^ 

For denominator values close to the value of 1 .0, an alternate equation can be used. 

(j) e = 4z^l + 4z 2 j 

Using the new equation for (j) e for circular arc camber, the corresponding lift and pitching- 
moment coefficient equations are: 

C / =4ttz(i + 4z 2 ) 
c mcl4 =-^(l + 4s 2 ) 

The factor |l + 4z~ j indicates that for 2% camber, the increase in the aerodynamic 
coefficients is of the order of one-sixth of one percent of the original value 
|l + 4z 2 j = 1 .0016 . This small amount is also the difference in the aerodynamic results 
between airfoils using circular arc and parabolic arc camber. 
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Pitch axis and moment center location Pitch axis and moment center location 

Figure 1. Thin airfoil theory results 



Figure 2. Effect on pitch damping of the number of vortices representing the airfoil 
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Pitch axis and moment center location 
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2L , chordwise position 

Figure 3. Comparison of discreet c 

vortex calculated values with thin Figure 4. Chordwise pressure distribution 

airfoil theory results fora = l/2ir 



0 0.2 0.4 0.6 0.8 1.0 


.£. , chordwise position 

Figure 5. Comparison of chordwise pressure 
distributions for two methods (a = l/2jt rad) 



-g-, chordwise position 


Figure 6. Chordwise pressure 
distribution for qc/2V = 1/it 
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a) Airfoil represented by two vortices 
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Figure 7. Sketches used in airfoil drag computation 
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Figure 8. Sketches illustrating locations of the incremental vortices for airfoil and 
wake with increasing time (Single resultant vortex represents the airfoil) 
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Figure 9. Sketches illustrating an improved arrangement (Two resultant vortices 
represent the airfoil) 
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Figure 10. Lift indicial response due to angle of attack 
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Figure 11. Pitching moment indicial response due to angle of attack 
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a) Effect of N 
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Figure 12. Response of induced drag function due to a step input in <x 
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Figure 13. Lift derivatives for sinusoidal oscillations in angle of attack 
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Figure 16. Moment - displacement trace for a sinusoidal plunging oscillation 
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c. Induced drag displacement trace for sinusoidal plunging oscillation 


Figure 17. Development curves for induced drag information 
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Figure 19. Pitching moment derivatives for sinusoidal oscillations in pitch rate 
holding angle of attack constant at zero 



c/2 

Figure 20. Induced drag function step input response to pitch rate for various pitch axis 
locations with N=4 (angle of attack held constant at zero) 
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Figure 22. Calculated aerodynamic derivatives for pitching moment of a thin 
airfoil oscillating sinusoidally in pitch about an axis on the chord 
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Figure 23. Comparison of results using Wagner's equation with 4-vortex 
model results 
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Figure 24. Comparison of Theodorsen prediction with 4- vortex model results 
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Figure 25. Comparison of TR 1108 results (pitch axis 0.37c) with experimental 
results corrected by lift curve slope ratio 
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Figure 26. Comparison of TR 1108 results (pitch axis 0.37c) with 4-vortex 
model predictions 
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Figure 27. Comparison of TR 1108 combination results (pitch axis 0.37c) with 4- 
vortex model predictions 
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Figure 28. Airfoil steady state lift and moment derivatives due to flap deflection for 
various values of the ratio of flap chord to airfoil chord (plain flaps with gap sealed) 
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b) Due to flap deflection 


H. 

c 


Figure 29. Flap lift to airfoil lift comparison for steady state conditions ( 
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Figure 30. Steady state flap hinge-moment derivatives as a function of the ratio 
of flap chord to airfoil chord (plain flap with gap sealed) 
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Figure 31. Lift derivative values due to a step input for the starting condition 

l?t vt 

(— = 0) and the terminal (steady-state) condition (— = oo ) for flap attitude (i.e., 
/2 /2 

flap angle of attack) 
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Figure 32. Moment derivative values due to a step input for the starting 

l?t Vt 

conditions (— =0) and the terminal (steady-state) condition (— = oo ) for flap 
/2 /2 

attitude (i.e., flap angle of attack) 






















Figure 37. Values of flap hinge-moment derivatives for the start (— = 0) and 

/ 2 

termination (— = oo ) conditions for step inputs in flap angle of attack and flap 

/2 

pitch velocity 
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Figure 38. Comparison of airfoil lift derivatives with the airfoil oscillating and 
with only the quarter-chord flap oscillating 


1c 1c 



Figure 39. Comparison of airfoil lift derivatives with the airfoil oscillating and 
with only the half-chord flap oscillating 
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Figure 40. Lift-displacement trace for an airfoil oscillating about the mid-chord 
with flap undeflected and in combination with a quarter-chord flap oscillating 
with the same amplitude and frequency 


50 



1.5 



a 

a 0 

Figure 41. Lift-displacement trace for an airfoil oscillating about the mid-chord 
with a quarter-chord flap oscillating with one-half the amplitude and twice the 
frequency of the airfoil motion 


51 




Figure 42. Lift-displacement trace for an airfoil oscillating about the mid-chord with 
a quarter-chord flap oscillating with the same amplitude and twice the frequency of 
the airfoil with 4>f = -60° 
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Figure 44a. Circular arc geometry 
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Figure 44b. Free stream flow normal component 
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Figure 44c. Camber geometry 
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Figure 45. Geometry for two circular arcs (max camber at 0.4c) 
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Figure 46. Geometry for two circular arcs (max camber at 0.6c) 





Figure 47. Lift and moment coefficients versus maximum camber location 
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Figure 48. Chordwise load distribution for three maximum camber locations 
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Figure 49. Camber geometry for NACA 4-digit series with max camber at 50% chord 
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